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Abstract 

We investigate superradiant cascade emissions from an atomic ensemble driven by two-color 
classical fields. The correlated pair of photons (signal and idler) is generated by adiabatically 
driving the system with large-detuned light fields via four-wave mixing. The signal photon from 
the upper transition of the diamond-type atomic levels is followed by the idler one which can be 
superradiant due to light-induced dipole-dipole interactions. We then calculate the cooperative 
Lamb shift (CLS) of the idler photon, which is a cumulative effect of interaction energy. We study 
its dependence on a cylindrical geometry, a conventional setup in cold atom experiments, and 
estimate the maximum CLS which can be significant and observable. Manipulating the CLS of 
cascade emissions enables frequency qubits that provide alternative robust elements in quantum 
network. 

Keywords: Superradiance; Cold atoms; Cascade atomic transitions; Two-photon spontaneous 
emissions; Quantum telecommunication 


1. Introduction 


Superradiance |[ll] is a coherent and collective radiation from a multi-atom system. The ra¬ 
diation intensity is proportional to the square of particle number along with a shorter time vari¬ 
ation, which conserves radiation energy. This collective radiation originates from a common 
light-matter interaction, through which spontaneously emitted photons can be reabsorbed by the 
atoms Hi if they are close to each other, and interatomic phase correlations build up 101 . This 
spontaneous emission decays in an enhanced rate proportional to the number of particles, which 
is very different from spontaneous emissions of independent and uncorrelated single atoms ||3l. 

Many theoretical works from different perspectives have investigated the superradiant emis¬ 
sion. The microscopic mechanism of superradiance is shown to be related to induced dipole- 
dipole interactions la [Zt] between two atoms, in which spontaneous decay rate and frequency 
shift fli depend on their spatial separations. The canonical review articles on the frequency 
shift including cooperative Lamb shift (CLS) and Lorentz-Lorenz shift from a local field correc¬ 
tion in the extended medium of a slab or a sphere can be found in the references of EEh. 
A master equation is formulated imB to study the statistical property of superradiance lll3ll . 
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and a diagonalization of coupling matrix QQ reveals the essential information of superra¬ 
diance decay constants. One of the characteristics in superradiance is its directionality of emis¬ 
sion il4ll6l] . and there is also a semiclassical treatment incorporating propagation effect, which 
indicates a threshold condition of cooperative emission in an extended medium OlTlIlSll . An al¬ 
ternative approach of quantum trajectory method ifl^ E(il] for superradiance tries to unravel the 
physics of successive sub- and superradiant photons, and positive-P phase-space method simu¬ 
lations show an enhanced decay in the second order correlation function by including quantum 
fluctuations in a diamond-type atomic system Olll . 

The preparation of single-photon absorption in the atomic ensemble followed by a collec¬ 
tive ll22n and directional emission ll23n raises the interests of single-photon superradiance. This 
singly-inverted system reduces Dicke’s full eigenstates of spin-1/2 system into only N state 


from a specihed initial state. The 


bases, and shows dynamical evolutions |24, 25, 2^ 27, 
initially symmetrical excited state for two-level atoms radiates superradiantly while other less 
symmetrical excited states radiate in a slower rate with smaller probabilities 11241] . Signihcant 
collective Lamb shift in a spherical geometry is also calculated in such singly-inverted system 
i^ . To correctly account for the cooperative effects in a dense medium, it is crucial to include 
the counter rotating-wav e-appro ximation (RWA) terms ll27ll rather than excluding them as in¬ 


complete treatments ll23L l24125L l26l] . Recently some experiments demonstrate superradiance or 


CLS in a variety of atomic systems, including atoms in aplanar cavity S, atomic vapor layers 
1 3in . cold atoms 13211 . and cascade atomic systems ll33ll34|] . 

The cascade atomic system provides a source for telecommunication bandwidth in its up¬ 
per transition 1^ . This correlated photon pair can realize long-distance quantum communica¬ 
tion e nabling a low-loss quantum repeater ll^ in the DLCZ (Duan-Lukin-Cirac-Zoller) 

protocols 11381] . In this paper we investigate the superradiance of correlated photon pair from a 
diamond-type atomic ensemble driven by two-color classical helds, and calculate the CLS in an 
optically-thick cylindrical atomic system. We solve for Schrodinger equations for such atomic 
system, and derive two-photon state function in section 2. Single atomic excitation adiabatically 
follows two excitation helds, and subsequently decays through cascade transitions via four-wave 
mixing (FWM). Counter-RWA terms in the Hamiltonian are required to correctly deduce the 
CLS which is a cumulative dipole-dipole interaction energy. In section 3 we calculate and study 
the CLS of two-photon state in a cylindrical geometry, a conventional setup in cold atom exper¬ 
iments. We conclude and summarize in section 4, and the inclusion of virtual photons exchange 
and adiabatic approximation of excitation process are detailed in Appendices A and B respec¬ 
tively. 


2. Two-photon state function 

2.1. Hamiltonian and Schrodinger equation of motion 

We consider an ensemble of N four-level atoms excited by two classical helds, and subse¬ 
quently the signal and idler photons are spontaneously emitted as shown in Fig. [1] Without loss 
of generality, we assume that these identical atoms distribute randomly in a cylindrical geometry 
with a uniform density as in conventional cold atom experiments 13311 . We use dipole approxi¬ 
mation of light-matter interaction and rotating wave approximation (RWA), and the Hamiltonian 
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Figure 1: (Color online) Four-level atomic system driven by two classical fields {Qa,b) spontaneously emitted photon 
pair (a^ ^). Single and two-photon detunings are Ai^ 2 - Schematic collective couplings between N state bases are denoted 
by arrows, and slightly shifted energy levels (not to scale) are due to light-induced dipole-dipole interaction energy. 


in interaction picture reads 

N N , 

Vi = -hA, ^ |1)^(1| - hA 2 Y. |2)^(2| - + h.c.} 

^ — 1 fl—l 

ks ^ -^s 

+ E 5k.(ek„A. • - h.c.}, (1) 

ki.Ai 

where h.c. denotes Hermitian conjugate. Note that the RWA Hamiltonian is valid for deriving 
state functions to calculate transition probabilities while non-RWA or counter-RWA terms are 
required to account for a complete derivation of CLS (see in Appendix A). The collective dipole 
operators are defined as 

PI ^ E PI ^ E 

fl fl 

si ^ E |2)/.(3|e*‘^--V ll ^ E (2) 

fl 

where single photon detuning Ai = LUa — wi, two-photon detuning A 2 = tUa + Wf, — u; 2 , and 
transition frequency W23 = u;2 — W3. Driving Rabi frequencies are Hq = (l||d||0)£^(fca)//i, fib = 
{2\\d\\l)£{kb)/h, and coupling coefficients are pk^, = (3||d||2)f (ks)/li, g^. = (0||(i||3)f (ki)//i. 
The double matrix element of the dipole moment d is independent of the hyperfine structure, and 
£{k) = ^/hkc/{2eoV) where V is the quantization volume. Polarizations of signal and idler 
fields are Ck^ ^.a^ i, and the unit directions of dipole operators are ds,i respectively. 

In the limit of large detuned and weak driving fields, Ai » VNlflal we consider only 
single excitation, and neglect spontaneous decay during excitations (Ai ^ Ti, A 2 ^ r 2 ). We 
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may express the state function as 


N N N 

IV’(O) “ ^(01^5 vac) + ^ ^ Il/j,j vac) + y] vac) + ^ lk,.A.) 

fl—1 ^ — 1 fl—l s 


+ X] lk,,A,, Iki.Ai), 


(3) 


where the indices s = (kg, Ag), i = (k^, Ai) for light fields, \m^) = for atomic 

levels m = 1, 2, 3, and |vac) is the vacuum photon state. These bare states make a complete 
basis of collective excitations for identical particles, which describe a complete cycle of a single 
atom following the excitation and spontaneously emitted fields. We apply Schrodinger equation 
= Vi{t)\ip(t)), and the coupled equations of motion are 




iAf, = - AiA^, 

iB^ = - *y]5k.(ek.,A, • 

S 

- * X 5k. (efc..A. • d*)e ^^■■■-, 

i 


(4) 

(5) 

(6) 

(7) 

( 8 ) 


In the next subsection, we proceed to solve for two-photon state function from the above 
equation of motion in the bare state basis. 


2.2. Two-photon state in bare state basis 

Firstly we solve for signal photon state function by substituting Eq. (O into Eq. (|7]i, and the 
time derivative of {t) becomes 


C^{t) = gliel ■ - X X l^d^kk^A, • d* 

ki,Ai 

X [ (9) 

Jo 


We use Weisskopf-Wigner approach which assumes a weak-coupling limit and implic¬ 
itly a Markov approximation W3 ^ Ea with an intrinsic relaxation time set by l/Ea (see below). 
In addition the secular approximation suggests an ensemble size (Lg) limit of Lg/c « l/Ea 
for the approach to be valid S0]- The summation of field modes in the above leads to the 


4 


frequency shift and decay rate which are 


Jo 


i: i: i9.i‘ki<.,». 

i' ki,Ai 

= — (t) [TrS{uJi - uj^) - iP.Y.{^uJi - ^ 3 )“^] , 


= J2^sit) 


-F^AO - 


(10) 


where ^ and = |r^ — rj^j. involves the Lamb shift (i/ = jS) and collective 

radiation shift from dipole-dipole interaction (p ^ v) which we describe in details in Appendix 
A. The other frequency shift of derived by including non-RWA terms as shown in 

Appendix A. 

Therefore the differential equation of probability amplitude of C^{t) becomes 


^ cAt)F,A0 +* E ‘^c:{t)G,AA + ^c^m 

V v^fj, 


( 11 ) 


where Cl = r 2 /(r 3 / 2 ), and fl is Lamb shift. Ff^^, and G^i, represent the collective contributions 
to spontaneous decay rate and frequency shift, which are induced from the couplings between N 
state bases as shown schematically in Fig. [T] 

The differential equation of Eq. (fTTT) can be expressed in a matrix form as an eigenvalue prob¬ 
lem, and we can solve for by a similarity transformations U = {ui,U 2 , ■■■, mat) comprised 

of eigenvectors Ui, 

l,m 




where Xi and Ui are eigenvalues and eigenvectors respectively of N x N matrix M, 


M = -- 


/ 1 

Fi 2 + i2Gi2 

V Fin + i2GiN 


Fi2 -f i2Gi2 
1 


Fin + i2GiN ^ 

1 / 


(13) 


Note that we have absorbed Lamb shift into the optical frequency W 3 , and the properties 
Fai 3 = Fpa, Gap = Gpa h^ve been used. The diagonal element demonstrates a single atomic 
spontaneous decay, and off-diagonal elements represent the couplings between atoms. The cross 
couplings are in similar order of magnitude indicating highly dynamical interactions in the sys¬ 
tem. 

We then substitute Eq. (fT2l i into Eq. (O, and derive the two-photon state function in the bare 
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state basis 


DsAt)=9:9UeLx.-l)K^,x,-d 


V7‘ 


-1^23-A2)t" 


^,l^m 


— zks 


(14) 


In the adiabatic approximation of the excitation process described in Appendix B, we may derive 
Bm{t), and the single and two-photon probability amplitudes become 


CAt) = 9 *A^L,x. ■ ds) f ^ 

l,m 

(15) 

DsAt) = • 1) f f 




(16) 


where the time evolution of two-photon state incorporates two-time integration of pulse shapes 
b{t) = /(4Ai A 2 ) which results from adiabatic driving process (See Appendix B). 

The cascade two-photon state is expressed of discrete sum of N eigenvalues A; shown in 
the time evolution and sum of eigenvectors sandwiched by phase factors induced by four fields 
interacting with the system. In the above treatment where we use the bare state basis, the signal 
spontaneous emission comes from the adiabatic transfer of detuned two-color driving fields, and 
the idler spontaneous emission comes from N intermediate excited states. These solutions of 
probability amplitudes demonstrate a dynamical coupling between N atoms, and two-photon 
state is generated from all possible excitations of one of the atomic ensemble. In the next section, 
we rotate the coupling matrix M to the new basis, and we find the superradiant idler photon with 
its decay constant coming from one of the eigenvalues. 

The superradiance (A; > Fa) and subradiance (A; < Fa) are embedded in these eigenvalues 
\i which depend on the density and geometry of atomic ensemble. For a low density where 
^ l/A, Fap ~ 0, Ga /3 ~ 0, and the coupling matrix M becomes an identity matrix times an 
overall constant —Fa/2. The transformation matrix U is then also an identity matrix. Therefore 
the two-photon state from N diamond-type atoms behaves as no difference from a single atom. 
On the contrary in Dicke’s limit where the dimensions of ensemble are smaller than radiation 
wavelength, Fajs ~ 1, and Gap becomes divergent when rap 0. Without considering this 
divergent CLS, one of the eigenvalues becomes Re(Ai) = — AFa/2 proportional to the number 
of particles, reminiscent of Dicke’s state |( = N/2,m = 1 — N/2) llllst] where m is the half 
difference of atomic populations in the excited and ground states. 

In the next section, we construct the basis in terms of the symmetrical and N-1 unsymmetrical 
states. In this basis we are able to investigate the superradiance and CLS of the two-photon state. 


2.3. Two-photon state in phased symmetrical state 


Here we investigate the collective radiation decay and energy shift of two-photon state in the 
basis constructed from the phased symmetrical state ill which provides a preference for the 
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system in the limit of large number of atoms. Consider a unitary transformation S where a phase 
factor in an extended ensemble is introduced into basis ll24ll . 


!</>,) = 




1 

7n 


flfi = + 


^/N 


i + i/Viv\ 
TV- 1 I 


’fiN 


'i + 1/Vn 

TV- 1 




'ill 


(17) 

( 18 ) 

Si^N, (19) 


where k = — kg, and the extra phase factor introduced above is from the observation in 

Eq. CSl). 

The phased symmetrical state in the new basis is 

= (20) 

and the other — 1 unsymmetrical states are 

\^i) = -^|3)Ar|0)f^^i (21) 

where I = 1,2,..., TV — 1, and they are orthogonal to each other and normalized to one. 

In this new rotational basis, the signal state function has the identity 

N N 

E E c’r(oi3^,ik..A.) = E E 

I-L —1 fcs,As / — I /Cs,Ag 

where these coefficients are related by the unitary transformation matrix that 

and Ci = 3*^0^. A new set of equation of motion in this basis can be derived, and we have 

for the equation of C*, 

Ci = "EE (22) 

I' 


where we have used the property of = x/NSkN, and M is the coupling matrix in bare 

state basis. The coupling matrix in this new rotational basis becomes (let + 


M = {SAS'<f 

_ "Ylii^i, f'ivfliiKiiviO 

2 : 

\ ThlL^U fNvflllKlwiC) 




(23) 
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where the properties = '/NSin, fivfi'u = Sw, M = are used, and T means 

transpose. The diagonal element An, i < N — I has the order of 1, which indicates a spon¬ 
taneously decayed single atom, and n involves the superradiant decay rate and CLS of the 
symmetrical state. 

Similarly we may diagonalize the above matrix with a similarity transformation matrix U, 
and the solutions for probability amplitudes of single and two-photon states become 


Ciit) = ■ 1) 


'U, 


1 ^ ’ 


(24) 


(25) 


M 1,3 


where Ak = -f kf, — kg — indicates FWM mismatch. The state functions show a specific 
dependence on Ath column of U~^ which signifies the difference of the new rotational basis 
from bare state one. 

One way to simplify the solutions is to examine the matrix elements of A. In the limit of 
large number of atoms. Aim and A^i with z ^ A are negligible in the order of compared 
to Aattv- It can be shown to use a Gaussian density distribution in continuous limit ll23l with 
approximation of fig Ri N~^ — 5i^ where I A ^■ Ths order of magnitude specifies the coupling 
strength between the phased symmetrical state and the other TV — 1 less symmetrical ones. One 
can show the cross-coupling elements of iV — 1 non-symmetrical states are not as negligible 
as Aifi or Ajsii. It is estimated as less than 1 /compared to A^n, but it does not 
come into play with our two-photon state due to the weak coupling constants of Ain and To 
lowest order of simplification of matrix A, we legitimately set iVth row and column of A to zero 
except for the element Ann and we find that Uin = Sin and U~j^ = Sjn which can be regarded 
as decoupling from this symmetrical state \(j>N)- The two-photon state therefore becomes 


t rt' 


DsAt) = 5r5:(4..A. ■ rfO(4,.A. • i) E / / 




dt”dt' -U>23-^2)t" 

(26) 


The above is one of the central results in this paper. The condition of FWM is embedded in 
the summation of four-wave phase factor gjAk r^ which is maximized when Ak = 0. The 
cascade emission has the most significant contribution of the phased symmetrical state where 
A N represents the timescale of spontaneous decay and CLS for the idler transition. Note that the 
cooperative effect for this correlated photon pair appears only on the idler transition, which is due 
to the large-detuned FWM driving process that the atomic system is adiabatically prepared to the 
upper singly-excited state. In deriving this two-photon state function with induced dipole-dipole 
interactions in the idler transition, the superradiant decay rate and CLS of the idler photon have 
no difference from the emission of a phased symmetric state in a two-level atomic system. 
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Figure 2: (Color online) (a) A superradiance decay factor + I’s dependence on height h and radius a in a cylindrical 
geometry. The enhanced superradiance factor has a stronger dependence on the height, which can be seen in the sectional 
plot of (b) while it saturates as radius increases as shown in (c). The legends in (b,c) denote the values of log(a/A) and 
log(h/A) respectively. D1 transition of the idler photon is chosen for Rb atoms, and the atomic density is 8 x 10^^ 
cm“^. 


The real part of Aat can be expressed in terms of geometrical constant [I16I] of a cylindrical 
ensemble. 


Rq{Xn) = —- 


1 ^ 




£3 

' 2 




-ik-(r^-r„) 


- -N / dOa — 

2 ./ 'Stt L 


1-Ch-df r(k,k 3 ) = -^(iV/i + l), 


(27) 


where in the second line we use = 1. fla is the solid angle of idler photon with wavevector 
ka, and r(k, ka) = ^ g*(k-k 3 ) (r;.-r„) yvJjJqJj ensemble average with randomly 


distributed atomic positions and r^. In deriving Eq. (l27l i we first substitute 
for / df 23 |f[l - (fca • d)2]e*k3.(r,-r.) ||^ as in Eq. (fTOl) where the property of polarizations 
^k.Afk.A = 1 — kk/|kp is used with kk in a dyadic form 1^ . Then we go to the continuous 
limit by replacing with y f cPr, and relate to /2 of a cylindrical geometry 10. We note 
that the excitation (absorbing two-color classical fields with an emitting signal photon) propa¬ 
gates along the height of the cylinder, which is k. The geometrical constant p for a cylindrical 
ensemble (of height h and radius a) therefore is 0 


d = 


da;(l 


6(V-1) r 

NA^H^ Ai (l-a;)2(l-a;2) 


sin 


H{l-x) 


’^1 


A^/T^: 


(28) 


where H = k^h and A = k^a are dimensionless length scales, and circular polarizations are 
considered lllhll . Ji is the Bessel function of the first kind. As shown in Eig.©, we demonstrate 
how superradiance decay depends on length and radius of a cylindrical ensemble. 

The real part retrieves Dicke’s result when atoms are confined within the radiating wave¬ 
length and = 1 — 1/iV. When the average separation of the ensemble is much larger than 
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Figure 3: (Color online) Cooperative Lamb shift (CLS) in a cylindrical geometry, (a) 3D plot of CLS on height h and 
radius a, and (b,c) cross sections on h and a respectively with the values of a and h indicated in the legends. The atomic 
density and other parameters are the same as in Fig|2] Short radius behavior involves a sharp increase which can be seen 
in Fig. |4lb). 


emission wavelength, dipole-dipole interaction becomes negligible, and we have single atomic 
spontaneous emission p. ^ 0. The frequency shift of this symmetrical state can be significant 
if atoms are close to each other. It is a cumulative contribution from dipole-dipole interactions 
which are oscillatory in space as shown in Appendix A. 


3. Cooperative Lamb shift of the phased symmetrical state 


In this section, we calculate the cooperative Lamb shift of phased symmetrical state, which 


IS 


Im(AAr) = - 


-ik-(r^-r„) 


dk 


—^PV 

Nkl ■ ■ y_oo 27r fc - fca 
' dk k^ 




-ik-(r^-r„) 


^3 i-oo 2TTk-k3 


N^i{k), 


(29) 


where P.V. denotes principal value, and we note that geometrical constant p{k) has a momentum 
dependence in the integral. Interestingly the CLS Im(AAr) relates to the integral of spontaneous 
decay rate proportional to Np{k), which is a Hilbert transform lld^ if we put back the Lamb shift 
into our expression. The Lamb shift is an intrinsic level shift induced by vacuum fluctuation while 
the CLS depends on the size and geometry of the atomic ensemble indicating a very different 
physical mechanism. 

To calculate the integral, we may start from the analytic form of fi{k), and re-express Eq. 


10 


























E 



0.03 



a (mm) 


0.05 
a (mm) 


Figure 4: (Color online) Collective Lamb shift of alkali metallic atoms as a dependence of radius a. The inset in (a) 
shows the crossing of red-blue shifts for different alkali metals, and short radius dependence is shown in (b). The height 
is chosen as = 3 mm. The atomic density and idler photon transition are the same as in Fig[2] 


( l29l l by renormalizing k ^ k^k. 


Im(A7v) = TaRV. f 

J —( 


dk k^ 6{N — 1) 
^k-1 k^A^m 


sin^ \}^kH{l — x)] J? \kA\/\ — 

X / l + ic^ - ^ ^ - Ux. (30) 

J-i (1- xy{i- x-^) 

The function of integral over x changes slowly compared to the change of k, and we may move 
it outside of the integral similar to the nonperturbative treatment in the radiation level shift 0. 
Taking k = 1 for the integral of x where the integral of k is most appreciable, we have 


Im(Ajv) = NfL{k3)r3 X RV. [ 

Jk. 


T 


1 1 
k — 1 k + 1 


(31) 


where we introduce the infrared (km) and ultraviolet (kM) energy cutoffs fflEl] to the integral, 
which represent the lowest and largest energy scales in the system. Then we have the CLS in the 
leading order of fcm /^3 and ks/kM, 


lm{XN)=2Nfl{h)Ts 


h 

h 


h 

kn 


(32) 


which indicates a redshift 0291 l3l[ 13211 for the level |3) if we let kM oo. In Fig. CLS shows 
a geometrical dependence on height and radius in a cylindrical atomic ensemble. We use km = 
j ^'kB?L as a long-wavelength momentum cutoff for the system, and use atomic radii 1421 
as a measure for maximal momentum cutoff which is reasonable for we can not probe atom’s 
internal structure in our quantum optical treatment. For a needle-like cylinder where H ^ A, we 
have a redshift of several kHz that can be observable in conventional experimental setup. A blue 
shift is also shown in Fig. Ofor a disk-like geometry (A < H) which is due to the counteracting 
term from system energy cutoff in Eq. (l32l i. 

In Fig. 0] we show the CLS as a dependence of radius, and demonstrate the crossing points 
of red-blue shifts for other alkali metals in the subplot (a). Due to different D1 transition wave¬ 
lengths of alkali metals, sodium atoms reach the zero frequency shift at a smaller radius compared 
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to rubidium and cesium ones. By scanning through the geometrical dependence and comparing 
with different atomic species, we may determine the system energy cutoff. The introduction of 
energy cutoff means lacking knowledge of the system momentum distribution as in the calcu¬ 
lation of Eq. (E)- If the experiment shows a redshift as radius varies, it means a unnecessary 
introduction of the short-wavelength energy cutoff, which suggests that CLS is only relevant to 
the macroscopic ensemble length scales {H and A in our case). 

Note that when H or A approaches the condition when ^ 1, CLS is vanishing as can be 
seen in the general form of Eq. ( l29b where Np.{k) 0. In Eig. |5b), we specihcally show a 
short radius dependence where there is a sharp increase of CLS. The maximum appears at about 
a ~ 20 /rm, which can be signihcant if we have an even higher atomic density. An estimated 
of frequency shift of MHz is in the reach in the D1 transition of Rb atoms for atomic density 
p ^ 10^^ cm“^. We note that similar red to blue shifts in CLS is demonstrated in an ellipsoid 
geometry ll^ . 


4. Conclusion 


We investigate the superradiance and cooperative Lamb shift (CLS) of correlated two-photon 
emissions in a diamond-type atomic ensemble. The photon parr is generated from the cascade 
transitions through four-wave mixing of two classical excitation helds. The superradiant idler 
photon is well described by one of the eigenvalues of coupling matrix constructed by the phased 
symmetrical state. This symmetrical state is signihcantly relevant when the number of atoms 
is large. We specihcally analyze the geometrical dependence of CLS on a cylindrical atomic 
ensemble. We hnd a redshift in a needle-like geometry, and suggest it could be blue-shifted in a 
disk-like one. We further demonstrate a comparison of the CLS for different alkali metals, and 
propose to determine the system energy cutoff in experiments. 

We can take advantage of manipulating system’s geometry and density to control the CLS, 
which enables robust frequency qubits in quantum network. In the same spirits of using dual¬ 
species matter and light-frequency qubits from a ®^Rb-®^Rb isotopic mixture ll^ . we may pre¬ 
pare atomic ensembles with two different geometries or densities that generate entangled photon- 


pair states (a 


I.ial.i + a 


t At 

s,2“i,2 


2 )|vac) assuming the same excitation processes, where we denote 


(1,2) for two different frequencies and (s,i) for signal and idler bosonic helds as,i respectively. 
We may construct multimode frequency qubits without acquiring multiple atomic species, and 
control the qubits via optical means of atom traps which can be dynamical and efficient. The mul¬ 
timode frequency qubits may offer a speedup of entanglem ent g eneration as in temporal modes 
of rare-earth-metal ions 14^ 3, spatial quantum registers |£7, ^], and multiplexing quantum 
repeaters | 
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Appendix A. Interaction energy 

In this section, we derive the self-interaction energy shift for single and many two-level 
atoms interacting with a single photon. We demonstrate that a correct formulation of the inter¬ 
action energy, including Lamb shift and CLS, is derived when non-RWA terms are included in 
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the Hamiltonian. These non-RWA terms indicate the energy nonconserving process. It involves 
virtual excitations which are unphysical but crucial in determining the frequency shift. 

In interaction picture for single-atom light-matter interactions, we express the Hamiltonian 
as 


k 

-f |0)(l|ake-*(‘^+“i)* -f (A.l) 

where |0) and |1) are ground and excited states respectively, and high frequency parts 
are non-RWA terms. If the atom is in the excited state, it may decay to the ground state with an 
emitted single photon, whereas the atom in the ground state may be excited also by emitting a 
photon from non-RWA terms. The state function can be expressed as 

1^) = C'o(f)lO) + Ci(f)|l) + ^Ci,k(f)|l, Ik) + ^ Co.k(t)|0, Ik). (A.2) 

k k 

The coupled equations from Schrodinger equation are 

k 

Cl = 

k 

Ci,k = 

Com = (A.3) 


Solving for the probability amplitudes of atomic energy levels, we have 

Co = - ^ f Co{t')dt', 

k 


(A.4) 


The summation of spontaneous emission modes in the above for the ground (sign -b) and the 
excited (sign —) states becomes 

J duj^ \tt5[uj ± wi) — iP.V.(a; ± , 

where the part of delta function simply represents the spontaneous decay rate, T = 

with the dipole moment d for the transition. RV. denotes principal value. Moreover the Lamb 

shift can be identified as the difference of two level shifts, 

J du}^ [P.V.(u; - - P.V.(w -b = H. 
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Now for many atoms interacting with a single photon including non-RWA terms, the Hamil¬ 
tonian becomes 

+ |0)^(l|ake*‘^''’''e-*(‘^+“i)‘ -f (A.5) 

where p denotes the atomic indices. The state function can be written as 

|^) = Co(t)|0)+^Cf^W|l/.,lk) + ^Cr(f)|U 

+ ^ Co.k(t)|0, Ik) + ^ ^ 1., Ik), (A.6) 

k ^,k 

where the last unphysical state (two atomic excitations with one photon present) is introduced to 
couple back to the singly-excited state |@]. The coupled equations of motion for the probability 
amplitudes are 


k,/i 

= -* E 5ke-*(‘"-‘"^)‘Co,k - * E E 5ke-*(^+‘"^)*e*-'^ 

k k 

Cf.k = -i5ke*(‘^+‘^^)*e-**^''''^Co, 

Co,k = 

M 

<^2.k = -*flke*^“+‘^^)*e-*‘^''''Cf, 


^2M + (m ^ i^) 


(A.7) 

(A.8) 

(A.9) 
(A. 10) 

(A. 11) 


where again the ground state energy level has a Lamb shift contribution similar to single atomic 
one except that an extra factor of N, the number of atoms, appears when substituting Eq. ( IA.9b 
into Eq. (lA^ . which shows its intrinsic, not collective, property by vacuum fluctuations. 

Carefully calculating the contribution from nonconserving energy states in the excited state 
by substituting Eq. (lA.llI) into Eq. (IA.8b . we have 


k 


^2.k + ^ v) 




(N - l)Cf (f') -f Y 


(A. 12) 


where we can combine the first term in the above with the Lamb shift contribution from the 
ground state (a factor of N), and we deduce the positive frequency part of Lamb shift H+q,. 
By substituting Eq. (I A. 10b into Eq. (IA.8b for the energy conserving term and picking out the 
atomic index p, we deduce the other part of Lamb shift Again the Lamb shift is derived as 
H = — n+Q, which is the same as single atom case. 
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The second term in Eq. (IA.12b is the dipole-dipole interaction energy from energy non¬ 
conserving terms. Along with the energy conserving terms by substituting Eq. (lA.lOl) into Eq. 
(IA.8b and picking out the atomic indices other than p, we deduce the collective Lamb shift and 
spontaneous decay rate as 


f Y, Cl (^')e 


k-(r^-r„) 


- E cr(i) 




(A. 13) 


where the property of frequency shift, is used, and ^ = Iklr^j^, r^i, = |r^ — t^\ 

with the transition wave vector |k|. Here we show that the treatment of Schrodinger’s equation is 
equivalent to Heisenberg’s picture that and Gct,p are defined as 101 


FaAO = A l-(rf-^a/3) 


sin^ 


1 - ?>[d-fapf 


= -{n-^ + 


1 - {d-fapf 


COS^ 


-f 


1 - 3{d-fap)" 


cos ^ sin ^ 

sin ^ cos ^ 






(A. 14) 

for a ^ P, 
(A. 15) 


where d is the unit direction of electric dipole. The dependence of in the above earns the 
name of dipole-dipole interaction which is induced by the common light-matter interaction and 
radiation reaction. 


Appendix B. Adiabatic approximation of excitation process 


In this section, we derive the probability amplitudes of the excitation process in the adiabatic 
approximation. Eirstly, we substitute the signal photon state of Eq. (fT2b into Eq. and the 
summation of signal field modes in Eq. (|6]l becomes 


E 

k.,A„ 


/ dt'e 


■'piGks-‘^23-^2){t'-t) 'y ^ giks-(r^-r„ 




l,m 

To 

= -f Lamb shift term, (B.l) 


where = (w 23 + A 2 )|rn — rm|/c. The summation of the signal field modes is determined by 
the fast oscillating exponential factor of optical frequency, which is valid for the eigenvalues A; 
<C a ;23 -f A 2 . In the last step of derivation, we may absorb Lamb shift into the signal transition 
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frequency, and the above result indicates that the upper excited state radiates in a single atomic 
decay rate. 

In the limit of large detunings, 

|Ai|, IA2I ^ fla, r2, 


we can solve the coupled equations of motion by adiabatically eliminating the intermediate and 
upper excited states in the excitation process. In Eqs. Q and we use integration by parts 
to express probability amplitudes in the first order of 1/Ai. Note that we allow time-varying 
excitation fields, and let A^(i:)= we have 


Mt) 



2Ai 


mw) , . 0 .1 ^ 

2Ai ^ 


1 

2 




(B.2) 


B^{t) = e*(^2-tir2/2)t 


2 

L2 




'2(A2+ir2/2) ^A2^’ 


(B.3) 


where the initial conditions B{—oo)=A{—oo) = 0 are applied. The adiabatic approximation 
requires the driving pulses to be smoothly turned on. Therefore in the first order of adiabatic 
approximation, we derive 


A^i{t) R 

2Ai 

(BA) 

£:(f)R 

^1, 

(B.5) 

Aw « 

'' 4 A 1 A 2 “ ^ ’’ 

(B.6) 


The above results show that the probability amplitudes develop by following the driving 
fields, and the ground state is approximately unity in the limit of large detunings. The AC Stark 
shift in the ground state can be ignored if N \fla{t')\^dt'<S^Ai which is also required for 

the validity of assumption of single excitation in our scheme. 
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